Spatiotemporal chaos and the dynamics of coupled Langmuir and ion-acoustic waves 

in plasmas 



S. Banerjee, 1 ' 2 '^ A. P. Misra, 3, [j] P. K. Shukla, 4 '[f] and L.Rondoni 1 

1 Department of Mathematics, Politecnico di Torino, Turin, Italy. 
2 Micro and Nanotechnology Division, Techfab s.r.L, Chivasso, Italy. 
3 Department of Physics, Umed University, SE-901 87 Umed, Sweden^ 
^Institut fur Theoretische Physik IV, Ruhr-Universitdt Bochum, D-44780 Bochum, Germany. 
(Received 11 January, 2010; Accepted 13 April, 2010) 

A simulation study is performed to investigate the dynamics of coupled Langmuir waves (LWs) 
and ion-acoustic waves (IAWs) in an unmagnetized plasma. The effects of dispersion due to charge 
separation and the density nonlinearity associated with the IAWs, are considered to modify the 
properties of Langmuir solitons, as well as to model the dynamics of relatively large amplitude wave 
envelopes. It is found that the Langmuir wave electric field, indeed, increases by the effect of 
ion- wave nonlinearity (IWN). Use of a low-dimensional model, based on three Fourier modes shows 
that a transition to temporal chaos is possible, when the length scale of the linearly excited modes 
is larger than that of the most unstable ones. The chaotic behaviors of the unstable modes are 
identified by the analysis of Lyapunov exponent spectra. The space-time evolution of the coupled 
LWs and IAWs shows that the IWN can cause the excitation of many unstable harmonic modes, 
and can lead to strong I AW emission. This occurs when the initial wave field is relatively large 
or the length scale of IAWs is larger than the soliton characteristic size. Numerical simulation 
also reveals that many solitary patterns can be excited and generated through the modulational 
instability (MI) of unstable harmonic modes. As time goes on, these solitons are seen to appear in 
the spatially partial coherence (SPC) state due to the free ion-acoustic radiation as well as in the 
state of spatiotemporal chaos (STC) due to collision and fusion in the stochastic motion. The latter 
results the redistribution of initial wave energy into a few modes with small length scales, which 
may lead to the onset of Langmuir turbulence in laboratory as well as space plasmas. 

PACS numbers: 52.35.Mw; 52.35.Fp; 52.35.Ra; 05.45.-a. 



I. INTRODUCTION 

The formation of envelope solitons, through the nonlin- 
ear interaction of high-frequency (hf) electric fields and 
low-frequency (If) ion- acoustic waves (IAWs) is one of 
the most interesting features, and is extensively studied 
problem in the context of turbulence in modern plasma 
physics. Formation of these structures not only plays a 
role in turbulence, but also for plasma heating, particle 
transport etc., and so their properties are of fundamen- 
tal interest. In plasmas, such envelope solitons are well- 
known Langmuir solitons that are hf Langmuir waves 
(LWs) trapped by the density troughs associated with 
If IAWs (see e.g., Refs. [THS]). 

During the past few years several attempts have been 
made to investigate the dynamics of solitons including 
dressed solitons [3], soliton collapse [7 , nucleation of 
cavitons [8], radiation of Langmuir solitons [2 , Landau 
damping in partially incoherent LWs [9] etc. An exper- 
imental observation for the excitation of IAWs has also 
been reported in a Langmuir turbulence regime by Prado 
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et al [TP] . Moreover, an increasing interest has also been 
found to study the Langmuir turbulence as a result of 
chaos [TTUT6] . When the electric field intensity is strong 
enough to reach the modified decay instability threshold, 
the interaction between LWs and IAWs results to 'weak 
turbulence', and then the LWs are scattered off IAWs un- 
der T e » Ti, where T( e ,i) is the electron (ion) temper- 
ature. On the other hand, when the field intensity is so 
strong that the modulational instability (MI) threshold 
is exceeded, the LWs are then essentially trapped by the 
density cavities associated with IAWs. Such phenomenon 
can occur frequently in plasmas. The interaction is then 
said to be in a 'strong turbulence' regime, in which trans- 
fer or redistribution of energy to higher harmonic modes 
with small wavelengths may be possible. In general, it is 
believed that some sort of chaotic process may be respon- 
sible for the energy transfer from large to small spatial 
scales jTT], [18] . This energy transfer may become faster 
when the chaotic process is well developed in a subsystem 
of the full set of equations [15] . 

One of the most important models in this context is 
the so-called Zakharov equations (ZEs) [T], which have 
been derived by assuming a quasineutrality as well as 
neglecting the ion-wave nonlinearity (IWN) due to den- 
sity fluctuation. However, the effects of dispersion due 
to charge separation (deviation from the quasineutrality) 
may be important in the radiation of IAWs, i.e., when the 
LWs resonantly interact with IAWs at the same speed 
[2j [3]. On the other hand, when wave amplitude be- 
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comes relatively large, the IWN can be important in the 
enhancement of ion-acoustic wave emission as well as in 
the excitation of many unstable harmonic modes at the 
initial stage of interaction, which later turns out to spa- 
tiotemporal chaos (STC). 

The basic purpose of this work is to investigate numer- 
ically the dynamics of coupled LWs and IAWs in presence 
of these important effects, and also to modify the results 
of Langmuir solitons not observed in the one-dimensional 
(ID) ZEs. We show that transfer of energy to fewer 
modes, indeed, occurs and it is faster when the chaotic 
process in a low-dimensional model of the full system is 
well developed by the IWN. Such low-dimensional model 
is constructed on the basis of three unstable harmonic 
modes [13] , and which depends on a particular regime of 
the wave number of excitation. Numerical simulation of 
this low-dimensional model reveals that transition from 
order to chaos is possible, when the largest length scale 
of the IAWs is greater than that of the soliton charac- 
teristic size. Solitons will then much be distorted by the 
I AW emission until a mostly chaotic state emerges. For 
details of the idea of constructing such model, readers 
are referred to the work in Ref. [T6] . 

Since the low-dimensional model is restricted to a few 
number of modes, and also it is valid for a certain vale of 
the wave number of excitation, one needs to investigate, 
e.g., numerically the the full space-time problem, espe- 
cially when this wave number is small enough to excite 
many unstable modes. We show that these solitary pat- 
terns, thus formed, ultimate lose their strengths after a 
long time through the random collision and fusion under 
the strong IAW emission. The STC state of the system is 
then said to emerge. However, the mechanism leading to 
such STC is still not clear. A few works in this direction 
can be found in the literature p~2j [16] [19] . 

The paper is organized as follows. In Sec. II, we 
present the governing equations, and analyze the evo- 
lution of Langmuir solitons in some particular cases of 
interest. A low-dimensional three- wave model is formed 
and its linear stability analysis is carried out in Sec. III. 
Numerical results for the temporal evolution of the low- 
dimensional model and the spatiotemporal evolution of 
the full system are presented in Sec. IV and Sec. V re- 
spectively. Finally, Sec. VI contains some discussion and 
conclusion. 



II. GOVERNING EQUATIONS AND THE 
LANGMUIR SOLITON 



The nonlinear interaction of hf LWs and If IAWs is 
governed by the following coupled set of equations [2j [3] 
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where E(x, t) is the slowly varying Langmuir wave elec- 
tric field normalized by y/l§irnok,BT e , with no denoting 
the equilibrium numbe r density, k B is the Boltzmann 
constant, e = ^m e /rrii is the square root of the electron- 
ion mass ratio and v(x,t) = n e i/no = n e — 1 is the 
plasma relative density fluctuation. The time t is nor- 
malized by the ion plasma period uj~ t 1 = 1/ ^ ^utiq e 2 /rrii 
a nd the space x by the electron Debye radius A^ e = 
^/k B T e /47rn e 2 . In Eq. ^ the terms with fourth-order 
mixed derivatives describe the effects due to charge sep- 
aration (i.e., deviation from the quasineutrality) . The 
third and fourth terms in the right-hand side of Eq. Q 
arise due to the expansion of log(l + v) up to the second 
order of v in the equation for If perturbation of electrons. 
This enforces us to assume that the relative density fluc- 
tuation associated with IAWs is not so small. 

We note that the analytic treatment of the Eqs. ([I]) 
and ([2| without IWN has been studied in detail for 
the modification of Langmuir solitons by Karpman et al 
[2j [3]. However, when the IWN is included in the sys- 
tem, the soliton properties may be changed. To look 
for this effect, we simply disregard the mixed derivatives, 
and solve numerically the system of Eqs. and ([2| by 
Runge-Kutta scheme. In the scheme, we approximate the 
spatial derivatives by the centered second-order differ- 
ence formulas, and use v = E = d 2 E/dx 2 = d 2 vjdx 2 — 
at the boundaries. The results are presented in Fig.l. 
When the ion-acoustic radiation is weak, we see that 
around x = the ponderomotive force dominates over 
the IWN, giving rise an increased amplitude of the enve- 
lope (highly correlated with the density depletion) . How- 
ever, in the other regimes of x, IWN dominates and 
makes the solitons wider. Later, we will see that as time 
progresses, the IAW emission increases by the IWN, caus- 
ing solitons to lose their strengths. Some particular cases 
can be of interest in the context of Langmuir solitons. For 
example, under the small amplitude condition, i.e., disre- 
garding the third and fourth terms in the right-hand side 
of Eq. Q , and assuming the quasineutrality (disregard- 
ing the mixed fourth order derivatives), the dynamics of 
small amplitude IAWs are then governed by 
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This equation when coupled with Eq. ([!]) describes under 
suitable renormalizations the classical ZEs [1 , which has 
been extensively studied in the context of Langmuir wave 
turbulence (see e.g., EB [HUH). 

Again, assuming a slow space-time response to the 
perturbations as x — >• e 1 ^ 2 (x — t),t —> e 3 / 2 £, i.e., ion 
p erturbatio ns move with the ion-acoustic speed, c s = 
\/k B T e /mi, and asuming the dependent variables to 
vary~ e, the IAWs can be described by the following 
modified Korteweg-de Vries (mKdV) equation [20H22] 
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FIG. 1: (Color online) The profiles of the Langmuir wave elec- 
tric field (upper panel) and the associated density depletion 
(lower panel) are shown [solution of Eqs. (1) and (2)]. The 
solid and dashed (red) lines respectively correspond to the 
case with IWN and without IWN after t = 10(cj~ 1 ). 
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The equation ([!]) together with Eq. Q describe the dy- 
namics of slowly varying Langmuir wave electric field cou- 
pled to the slow response of the density perturbations. 
These coupled equations predict a class of solutions with 
a bell-shaped Langmuir wave electric field profile arrested 
into the density cavity [2TJ [22] . 

We now turn out to the global behaviors of the system. 
To this end, we first renormalize the original variables 
(similar to the ZEs) as x — >• ex,t — >• et,v — >• v/e 2 ,E — >• 
E/e, and obtain the condition for MI as given below. 
Thus, Eqs. and Q reduce to 
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Clearly, the term proportional to e 4 (~ 10 -7 ) becomes 
smaller and can be neglected. However, the terms pro- 
portional to e 2 can no longer be neglected when the soli- 
ton attenuation becomes relatively larger. 



Next, the linear stability analysis of the perturbations 
oc ex.p(ikx — iuot) [where k (u) is the wave number (fre- 
quency) of modulation] for Eqs. ([5| and (|6| gives the 
following growth rate of MI. 
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where a = 1 - 



-k 2 e 2 , A : 



1 + 3/c 2 (3a/2 - e 2 |£ | 2 ) A and 



the MI sets in for k satisfying < k < k c 



III. LOW-DIMENSIONAL (THREE- WAVE) 
MODEL 

We note that the system of Eqs. ([5| and ([6| is, in 
general, multi-dimensional. However, there may be the 
situation in which a few number of modes is more active 
than the remaining modes. Such cases are quite typical 
not only for Langmuir wave turbulence, but are rather 
general for parametric instabilities of hf waves interact- 
ing with If ones close to the instability threshold. In 
this way, Galerkin expansions and truncations to a few 
normal modes are commonly used to describe the ba- 
sic features of the full dynamics by a low-dimensional 
model. However, we will see that the specific details of 
such model are highly dependent on the range for the 
basic wave number of modulation k. So, considering the 
dynamics of few coupled waves, we expand the envelope 
E(x, t) and the density v(x, t) as [T6] 
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where M = [& _1 ] represents the number of modes, p rn = 
p-rru Qui = n rn = n~m- Since in our case, M = 2, 
pi = 774 = n±i etc., we can express the fields as (see 
for details [IU [T3[15]) 

E = E + Ei cos(kx) = -/N sin {^j exp(i0 o ) 

+ V2N cos (|) exp(i(9i) cos (foe) , (10) 



v = N + u\ cos(kx) , 
where N (= n ) |^-i| 2 + 



(11) 

| Ei | 2 is the conserved 



plasmon number. Also, a, #o, #i and ri\ are each a func- 
tion of time t. Inserting Eqs. (10) and (11) into the 
Eqs. Q and and following, e.g., Refs. [Til ITTJl [T5] 
we obtain the following temporal system 
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- ^3k 2 — \pln\ cos if cot a^j , 



where X = [#1,2:2, #3, #4] is a column vector and Df(Q) 
is the Jacobian matrix at Q, i.e., 



hi = Uin\ + f/ 2 ni + Z7 3 , 
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where the 'dot' represents differentiation with respect to 
t, ip = 9o — Oi, and Ui, U2 and U3 are given by 

.. , > „ COS If , . o 2 \ 

t/i(a, = Gi sin 09 + cos a), 

sin a v 7 

E/2 (a, <p) = C2 + (C3 — Ui sin a) cos a, 
Us (a, y>) = C4 sin a cos ip. 

Here the constants Ci, C4 are 

iVfc 2 e 2 „ k 2 (f - JV) 
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The system of Eqs. dl2j)-(14) can be recast as a set of 



four ordinary differential equations (ODEs): 
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V 2e 
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±3 = £4 , 

±4 = ^1^3 + ^2^3 + % , 
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where J7i, C/2, C/3 are now functions of xi, #2, and we re- 
define the variables as a = #1, = #2, ^1 = ^3, ^1 = #4. 
The basic properties of the system can now be in order. 
It can be shown that the adiabatic limit of Eqs. (17)- 



(20) gives an integrable system, and hence nonchaotic 



[13J [15l [16] . Next, the local stability of the system about 
the fixed points (#10, £20, £30, 0) given by 



sinx2o = 0, cot #io = ± 



3k 2 



(21) 



±3C 3 k 2 ± V2C 4 + CW9/C 4 + 2x2. 



9Cife 4 
x/2 



±1 



9FT24 



^30 



0, 



(22) 



30 



where £30 ^ 0,cosx 2 o = ±1, can be examined by lin- 
earizing the Eqs. ([17])- (20), and solving the correspond- 



ing eigenvalue problem. The linearized set of equations 
in a neighborhood of Q = (xio, £20, x 30? 0) i s 
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The eigenvalue equation \Df(Q) — XI\ = then reads 
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± Ci cos 2 X10 (1 — cos #10 cot £10) 
_l ( n X30 

± Gi COSX10 ; 

V sin #10, 

x (sin(2#io) + cot #10 cos 2 #10) , 



A 2 = ±2x30 cos #10 , A 3 = A 2 cot xio . 



(27) 



(28) 



We numerically examine a set of fixed points satisfying 
Eqs. (21) and ([22|, and the linear stability of the system 



by the Routh-Hurwitz (RH) criteria. To this end, we use 
Carden's method to look for a real solution of Eq. (22) 
when expressed it as a cubic in a 



V2i 



30 



9/c 4 . The 

range of k and N for such real a, and hence for X30 are 
< k < 0.5 < y/2/3\E \ and 0.6 < N < 8. Thus, for a 
given set of values of k and TV, one can obtain a set of 
fixed points from Eqs. (21) and (22). Having obtained 



a set of fixed points, one can then numerically examine 
the coefficients of the depressed quartic [Eq. (25)] for 
< k < 0.5 and 0.6 < N < 8. We find that in these 
regimes, coefficients of A 2 and A are always positive, and 
the constant term changes sign as in Fig. 2. Here the 
shaded (white) region indicates the positive (negative) 
values of the constant term. The constant term is also 
positive for < k < 0.49 and 0.6 < N < 8. Thus, by 
the RH criteria Eq. (25) has at least one root with pos- 



itive real part when the constant term of the quartic is 
negative. The system of Eqs. (17)-(20) is then said to 



be linearly unstable. As an illustration, the eigenvalues 
for TV = 8, 0.498 < k < 0.5, are a pair of complex conju- 
gates, one real positive (> 1) and one negative real roots. 



5 




FIG. 2: (Color online) The plot of the constant coefficient 
in the eigenvalue equation (25) with the variations of N and 
k, showing its positive (shaded colored region) and negative 
(white region) values. The coefficient is also positive (not 
shown) for < k < 0.49 and 0.6 < N < 8. 



IV. NUMERICAL SIMULATION FOR A 
THREE- WAVE MODEL: THE TEMPORAL 
DYNAMICS 

From Eq. ([7]), the MI sets in only if k < k c = 
y^2/3\Eo\, and k = k c defines the curve along which 
pitchfork bifurcation takes place [23]. Also, the dynamics 
is really subsonic in k c /2 < k < k c where the MI growth 
rate is small [16]. As, k decreases from fc c /2, many more 
unstable modes (since M = with higher harmonics 

will be excited. In this case, even though the dynamics is 
not subsonic, a description with three modes should be 
relatively accurate for k < k c /2. Our basic assumption 
at this point is that the Langmuir modes for \m\ > 1 
are stable, and the ion-acoustic modes with \m\ =2 
remain as they are, already excited by the presence of 
| ra | = 1 Langmuir wave envelopes. This results a sys- 
tem of ODEs [Eqs. (17)-(20)] describing the temporal 
dynamics of the field variables. Notice, however, that al- 
though the Fourier modes are truncated, no assumption 
on the slowness of the dynamics is made. We would like 
to see how the temporal system behaves as k approaches 
k c from fc c /2, and the condition for the subsonic region 
is progressively relaxed by reducing k from k c /2 to some 
extent. Here k is not so small that the three- wave model 
becomes inappropriate, because smaller the /c, the larger 
is the number of unstable modes. In the latter case, full 
simulation of Eqs. ([5| and Q will be necessary to give 
more appropriate information about the interaction. 

Thus, for the temporal evolution, we numerically in- 
tegrate the system of Eqs. (17)-(20) by the Runge- 
Kutta scheme with step size h = 0.001 and initial values 
x\ = 0.1, £2 = 0.1, xs = 0.2, £4 = 0.3. By suitably vary- 



ing the system parameters k and N together with the ini- 



FIG. 3: (Color online) A periodic limit cycle in the low- 
dimensional model for k — 1.5, N — 6 and Eq — 2, showing 
that the three- wave dynamics is stable. 



tial choices, we have obtained different sets of parameters 
which can exhibit periodic as well as chaotic attractors. 
The latter are established by the Lyapunov analysis with 
at least one positive exponent indicating chaos. On the 
other hand, the regular or limit cycles may be described 
on the basis of their structures, amplitude oscillations 
with time or Lyapunov exponents with negative or zero 
values. 

For N = 6, E = 2 and k = 1.5 in k c /2 < k < k c , 
Fig. 3 shows a periodic limit cycle. This suggests that 
within the plane wave region, strong chaotic activity is 
completely absent, i.e., most of the Kolmogorov-Arnol'd- 
Moser (KAM) surfaces, especially the central fixed point 
remains unaffected by the nonlinear interaction. In fact, 
the absence of chaos is not only for these parameter val- 
ues, but also for some other values smaller or larger. It 
also explains that within this plane wave region no mas- 
sive redistribution of energy takes place. This informa- 
tion well agrees with the full simulation of the space- 
time problem [16]. As k is lowered from fc c /2, e.g., 
k = 0.5, Fig. 4 shows that the periodicity or integra- 
bility of the temporal system tends to become poorer 
exhibiting chaotic orbits. This can be verified by the 
set of corresponding Lyapunov exponents as shown in 
Fig. 5, where one positive exponent indicates chaos. We 
also find that as k becomes smaller, the largest Lyapunov 
tends to increase from a positive small value to a large 
one (> 1). So, it may be possible that as the basic wave 
number k is decreased further, strong chaotic features 
develop into the system. Thus, when k < k c /2, the sys- 
tem's periodicity breaks down. This is expected, since 
for smaller /c, the solitons get more affected by the ion- 
acoustic radiation until they are completely destroyed. 
In the strong chaotic process one can expect that soli- 
tons will deliver their energy to the remaining modes so 
as to gradually disappear after some time. However, 
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FIG. 4: (Color online) A chaotic phase portrait in the low- 
dimensional model for k = 0.5, N = 6 and Eo = 2, showing 
that the three- wave dynamics is unstable. 



FIG. 6: (Color online) Supercritical Hopf-bifurcation in the 
low-dimensional model for k — 0.45, N — 7 and Eo — 2, 
showing the bifurcation of an attracting periodic orbit. 
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FIG. 5: (Color online) The four Lyapunov exponents with 
respect to time (te) corresponding to the same parameter val- 
ues as in Fig. 4. One of the exponents is clearly positive, 
indicating that the three- wave dynamics is chaotic. 



one can no longer rely on the low-dimensional model for 
some smaller values of /c, as it gives rise an excitation of 
large number of unstable harmonic modes. This is the 
situation where full simulation of Eqs. ([5| and (|6| will 
give better understanding for the interaction. 

In addition to the above results, one interesting can 
also be to study the existence of Hopf-bifurcation in the 
Eqs. (17)-(20) by choosing k as the bifurcation parame- 



ter. Generally, Hopf-bifurcation is defined as the change 
in qualitative behavior of the system when a pair of com- 
plex conjugate eigenvalues passes through the imaginary 
axis. A pair of such roots can be obtained by numerically 
examining the eigenvalue equation (25). Thus, one can 



perform a series of analysis with different values of k. Fig- 
ure 6 shows the supercritical diagram corresponding to 
N = 7 and k = 0.45 in which an attracting periodic orbit 
bifurcates. The detailed analysis of such bifurcation is, 
however, beyond the scope of the present work. 



V. SIMULATION RESULTS FOR THE 
SPATIOTEMPORAL EVOLUTION 

In the previous section we have observed that as the 
value of k is lowered from fc c /2, the system exhibits strong 
chaotic features. So, in order to investigate the global 
behaviors of the solitons, namely the soliton interaction 
with strong IAW emission, formation of incoherent pat- 
terns due to collision and fusion, it is reasonable to con- 
sider the system of Eqs. ([5| and (|6| without the mixed 
derivatives (valid for long- wavelengths) . In our numer- 
ical scheme, we assume the spatial periodicity with the 
simulation box length L x such that L x = 2ir/k, the res- 
onant wave length. The integer M = representing 
the number of modes, is considered as a number of grid 
points, e.g., 256,512, 1024 and 2048 as the case may be. 
We choose the initial conditions as |T2l H9l 



n(x, 0) 



E [1 + bcos(kx)/L] , 
—2bEo cos(kx) / L, 



(29) 



where Eq is the amplitude of the pump Langmuir wave, 
b is a suitable constant, and L is of the order of 10 -3 
to emphasize that the perturbation is relatively small. 
The spatiotemporal system was advanced in time using 
the Runge-Kutta scheme with a time step dt = 0.001. 
We consider the grid size M = 2048 for lower values of 
/c, and M = 1024 for relatively large /c, so that x = 
corresponds to the grid position 1024 and 512 respec- 



7 



tively. In the simulation, we consider Eq — 2, b = 1 and 
e = y/l/1840, unless otherwise mentioned. 

Figures 7 (in presence of IWN) and 8 (in absence of 
IWN) show the profiles of the electric filed (upper panel) 
and density fluctuation (lower panel) associated with the 
IAWs at the end of simulation for k = 0.22, i.e., for the 
of M = = 4 harmonic modes. At the central part, 

we observe an excited electric field, \E\ ~ 9.5e highly 
correlated with the density depletion, v = — 81. 2e 2 in 
presence of IWN. Whereas in absence of IWN, these are 
\E\ ~ 9.375e and v = -84.4e 2 . Clearly, the effect of IWN 
is to enhance the magnitude of the wave electric field or 
to decrease that of the density fluctuation. From Fig. 9 
(contour plot corresponding to Fig. 7), we observe that 
the solitary master patterns are first formed by the mas- 
ter mode, and then begin to appear from the unstable 
harmonic modes. During the pattern formation, ions are 
usually driven out (by the resultant of ponderomotive 
force and the IWN) from the regions of the master mode 
as well as from some harmonic patterns, and then they 
form density distributions with v > (see the lower pan- 
els of Figs. 7 and 8) along the spatial axis. These ion 
density humps move stochastically on either side of the 
patterns, and do not arrest the Langmuir wave fields. 
We observe from Figs. 9 and 10 (contour plot corre- 
sponding to Fig. 8) that the pattern selection leads to 
three (< M = 4) harmonic modes initially peaked at 
x = L x /4e,L x /2e and 3L^/4e. These three modes are 
basically generated by the master mode and the unsta- 
ble harmonic modes. The system is then in the coexis- 
tence of temporal chaos (TC) and spatio-partial coher- 
ence (SPC) state similar to the case of plane wave region 
(k c /2 < k < k c ) explained in the previous section. 

However, as the IWN is relaxed, only two solitary 
patterns initially peaked at x = L x /Ae and 3L X /Ae are 
seen to collide with the master mode at t « 120/e and 
t « 170/e respectively (see Fig. 10). Also, the collision 
is so weak that it does not influence the soliton dynam- 
ics. The reason is that in absence of IWN and for rel- 
atively small Eo, the ponderomotive force was the only 
force to cause the weak ion-acoustic wave emission. As a 
result, though the collisions take place but are not suffi- 
ciently strong to fuse into a new one. On the other hand, 
when Eq is relatively large, e.g., Eq = 4, i.e., when the 
amplitude of both the initial profiles increase, one may 
observe (see Fig. 11) more collision and fusion due to 
strong I AW emission by the IWN. In this case, \E\ in- 
creases, but \v\ decreases as compared to the case of no 
IWN (Fig. 12). Figure 11 also shows that two pairs 
of solitary patterns initially peaked at the either side of 
x « 256/e, collide and fuse into two new patterns, one 
of which again collides with the other two master pat- 
terns peaked at x ~ 256/e and x ~ 465/e at t ~ 115/e 
and 165/e. After these collisions, there remain only three 
distorted patterns. The system is in TC, but still in the 
coexistence of TC and STC, since a few collisions will not 
suffice to cause STC [19] . 

Now, if k is further lowered, e.g., k = 0.065, many 
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FIG. 7: (Color online) The profiles of the wave electric field 
(| -Ei | = \E\/e, upper panel) and the associated density fluc- 
tuation (y\ = is/e 2 , lower panel) with respect to the space 
(X = xe) after time t — 200/e in the numerical simulation of 
Eqs. (5) and (6) for k = 0.22 and E = 2 with IWN. 



10 




-100' — 1 1 1 1 ' 

— 1 -5 5 10 

X 



FIG. 8: (Color online) The same as in Fig. 7, but without 
IWN. 
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FIG. 9: (Color online) Contours of \E(x,t)\ =const. with 
respect to the space (xe) and time (£e) corresponding to the 
evolution as in Fig. 7. This shows that the pattern selec- 
tion leads to three harmonic patterns. The system is in the 
coexistence of TC and SPC state. 
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13 that the first collision occurs at t « 57/e between 
two solitary patterns. These two patterns then fuse to 
form a new one, which again collides with other three 
master patterns initially peaked at x « 190/e, 325/e and 
700/e at t « 95/e, 120/e and 150/e respectively. There 
are also other collision and fusion that can be explained 
similarly. In general, as time progresses, the patterns 
seem to be much distorted in losing their strengths. The 
original 17 (Fig. 13) or 15 (Fig. 14) solitary patterns 
are then finally fused into a few incoherent ones. As soli- 
tons gradually vanish, irregular radiation increases due 
to stochastic motion. That is, as soon as a chaotic sub- 
set with less degrees of freedom is formed, the subset will 
then act as a pump delivering a net amount of energy to 
its neighboring modes. This energy transfer takes place 
because of the very random nature of the couplings. The 
system is then said to be in the STC state. In this case, 
for t > 100/e the spatial correlation function will ap- 
proach to a zero value [19] . A certain amount of energy, 
which was initially distributed among the 15 or 17 soli- 
tary waves, will now be transferred to a few incoherent 
patterns as well as to some stable higher harmonic modes 
with short- wavelengths. 

So, if initially there exist many modulational unstable 
modes with different modulation lengths to form solitary 
envelopes, collision and fusion among most of them can 
lead to the state of STC. There must exist a critical value 
of this wavelength or wave number k at which the tran- 
sition from TC to STC takes place. This value of k is 
roughly in 0.03 < k < 0.1. Comparing the pattern for- 
mation, collision and fusion in Figs. 11, 12 and 13, 14 
one finds that IAW emission is more stronger in pres- 
ence of IWN. Thus, in order to observe the interactions 
between LWs and IAWs at longer time scale and for long- 
wavelength limits, the effect of IWN is more important, 
which can no longer be neglected especially for relatively 
large amplitude waves. In this way, the term contain- 
ing the IWN in Eq. ^ would act as a correction to the 
dynamics of arbitrary amplitude coupled Langmuir and 
IAWs. 



FIG. 10: (Color online) Contours of \E(x,t)\ =const. with 
respect to the space (xe) and time (te) corresponding to the 
evolution as in Fig. 8. This shows that the pattern selection 
leads to three harmonic patterns. Two collisions occur at 
t ~ 120/e and t w 170/e. The system is still in the coexistence 
of TC and SPC state. 



solitary patterns are formed from the master mode and 
unstable harmonic modes (see Fig. 13 with IWN and 
Fig. 14 without IWN) by means of pattern selection. We 
observe that when IWN influences the dynamics, many 
unstable harmonic modes are excited compared to the 
other case. As a result, not only two solitary patterns 
collide (as mostly in the case of no IWN) to fuse into a 
new pattern, but more than two also collide repeatedly 
to form new ones which again collide with other new pat- 
tern^) fused by the other collisions. Notice from Figure 



VI. DISCUSSION AND CONCLUSION 

From the analysis of the previous sections, we find 
that the three-wave model can be a good approxima- 
tion [16 for the interaction of coupled LWs and IAWs 
in the plane wave region k c /2 < k < fc c , where the sys- 
tem exhibits stable oscillations. It can relatively be ac- 
curate in the region 0.3 < k < k c /2 where the system 
exhibits temporal chaos for a given value of the pump 
electric field Eq. For values of k < 0.3, which indicates 
the excitation of many (> 3) unstable harmonic modes, 
the low-dimensional model fails to describe the dynam- 
ics of the coupled waves. We have found that the three- 
wave model exhibits periodic orbits in 1.2 < k < /c c ,and 
chaos in 0.3 < k < 0.5 < k c /2 for a given initial pump 
Eq = 2 and a given constant N = 6. The space-time evo- 
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FIG. 11: (Color online) The same as in Fig. 9 for k = 0.22, 
but for relatively large Eo, i.e., Eo = 4, indicating that more 
collision and fusion take place due to IWN. 
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FIG. 12: (Color online) The same as in Fig. 11, but in absence 
of IWN, indicating that relatively a less number of collision 
and fusion take place compared to Fig. 11. 



FIG. 13: (Color online) Contours of \E(x,t)\ =const. with 
respect to the space (xe) and time (te) for k = 0.065 and 
Eo = 2 in presence of IWN. The pattern selection shows that 
seventeen solitary patterns, which were formed initially, col- 
lide and get fused into other incoherent patterns after some 
time. The collision is random and not confined between two 
patterns. The IAW emission also occurs, and the STC state 
emerges. 




lution, on the other hand, reveals that the system is in 
the state of TC for relatively a lower value of fc, i.e., in 
the regime 0.1 < k < 0.22, and the STC state emerges 
for 0.03 < k < 0.1. 

To conclude, we have investigated the nonlinear in- 
teraction of coupled LWs and IAWs by considering the 
effects of dispersion due to charge separation, and the 
nonlinearity associated with the ion density fluctuations. 
The latter can cause the excitation of more unstable har- 
monic modes, and can lead to strong IAW emission in the 
space-time evolution. The low-dimensional model shows 
that a transition from periodicity to chaos is possible for 
a suitable choice of the wave number of modulation. It 
basically predicts some basic features of the system as 



FIG. 14: (Color online) The same as in Fig. 13, but in 
absence of IWN. Notice that fifteen solitary patterns were 
formed initially. After some time they collide and fuse to 
form several new incoherent patterns. Here the collision is 
mostly between two patterns compared to Fig. 13. The STC 
state is also said to emerge in this case. 



well as suggests an approximate region of k for the exis- 
tence of regular and chaotic structures. The space-time 
evolution indicates that solitary patterns are first formed 
by the unstable harmonic modes. These modes are then 
excited by a master mode due to initial MI. As solitons 
experience oscillatory motions, they emit radiation due 
to irregular interaction with the ion-acoustic fields. Such 



10 



radiations are broad-band in nature, causing a growing 
number of modes to be involved in the stochastic dy- 
namics. Depending on the modulation length scale, the 
system may either be in the TC and SPC states or in TC 
and STC states. If few harmonic patterns coexist with 
the master pattern, and IAW emission remains weak, the 
system is still in the state of TC and SPC. However, if 
many solitary patterns are formed at the early stage due 
to MI, the system may experience the state of SPC. 
Over longer time scales, when solitons are affected by the 
free ion- acoustic radiation, collision and fusion among the 
patterns occur frequently as well as repeatedly to form 
several new incoherent patterns. The STC state is then 



said to emerge in which energy can flow faster the smaller 
is the wave number of modulation k. Since our model is 
in one spatial dimension, energy transfer is not a result of 
soliton collapse, rather it is purely generated by the non- 
integrable nature as well as chaotic aspects of the system 
[T6] . The results may have physical significance to related 
problems, such as turbulence in pulsar radiation where 
strong magnetic field plays an important role [24] . 
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